A comparison of spectral element and finite difference methods using statically refined nonconforming grids for the MHD island coalescence instability problem
نویسندگان
چکیده
A recently developed spectral-element adaptive refinement incompressible magnetohydrodynamic (MHD) code [Rosenberg, Fournier, Fischer, Pouquet, J. Comp. Phys. 215, 59-80 (2006)] is applied to simulate the problem of MHD island coalescence instability (MICI) in two dimensions. MICI is a fundamental MHD process that can produce sharp current layers and subsequent reconnection and heating in a high-Lundquist number plasma such as the solar corona [Ng and Bhattacharjee, Phys. Plasmas, 5, 4028 (1998)]. Due to the formation of thin current layers, it is highly desirable to use adaptively or statically refined grids to resolve them, and to maintain accuracy at the same time. The output of the spectral-element static adaptive refinement simulations are compared with simulations using a finite difference method on the same refinement grids, and both methods are compared to pseudo-spectral simulations with uniform grids as baselines. It is shown that with the statically refined grids roughly scaling linearly with effective resolution, spectral element runs can maintain accuracy significantly higher than that of the finite difference runs, in some cases achieving close to full spectral accuracy.
منابع مشابه
Adaptive mesh refinement with spectral accuracy for magnetohydrodynamics in two space dimensions
We examine the effect of accuracy of high-order spectral element methods, with or without adaptive mesh refinement (AMR), in the context of a classical configuration of magnetic reconnection in two space dimensions, the so-called Orszag-Tang vortex made up of a magnetic X-point centered on a stagnation point of the velocity. A recently developed spectral-element adaptive refinement incompressib...
متن کاملA Hybridized Crouziex-Raviart Nonconforming Finite Element and Discontinuous Galerkin Method for a Two-Phase Flow in the Porous Media
In this study, we present a numerical solution for the two-phase incompressible flow in the porous media under isothermal condition using a hybrid of the linear lower-order nonconforming finite element and the interior penalty discontinuous Galerkin (DG) method. This hybridization is developed for the first time in the two-phase modeling and considered as the main novelty of this research.The p...
متن کاملSpectrally formulated finite element for vibration analysis of an Euler-Bernoulli beam on Pasternak foundation
In this article, vibration analysis of an Euler-Bernoulli beam resting on a Pasternak-type foundation is studied. The governing equation is solved by using a spectral finite element model (SFEM). The solution involves calculating wave and time responses of the beam. The Fast Fourier Transform function is used for temporal discretization of the governing partial differential equation into a se...
متن کاملFree and Forced Transverse Vibration Analysis of Moderately Thick Orthotropic Plates Using Spectral Finite Element Method
In the present study, a spectral finite element method is developed for free and forced transverse vibration of Levy-type moderately thick rectangular orthotropic plates based on first-order shear deformation theory. Levy solution assumption was used to convert the two-dimensional problem into a one-dimensional problem. In the first step, the governing out-of-plane differential equations are tr...
متن کاملAdaptive Unstructured Grid Generation Scheme for Solution of the Heat Equation
An adaptive unstructured grid generation scheme is introduced to use finite volume (FV) and finite element (FE) formulation to solve the heat equation with singular boundary conditions. Regular grids could not acheive accurate solution to this problem. The grid generation scheme uses an optimal time complexity frontal method for the automatic generation and delaunay triangulation of the grid po...
متن کامل